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基于 动态 EM-SHSMM 的 异常 数据 下 设备 健康 预测 研究 
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(上 海 理 工大 学 管理 学 院 ， 上 海 200093) 


摘 要 : 针对 设备 退化 过 程 中 异常 数据 下 的 剩余 有 效 寿 命 预测 问题 ， 提 出 了 一 种 基于 动态 的 期 望 最 大 化 算法 (EM)- 分 段 
隐 半 马尔 可 夫 模 型 (SHSMM) 预 测 方法 。 首 先 ， 基 于 SHSMM 的 理论 框架 ， 采 用 期 望 最 大 化 参数 自 适应 估计 算法 估计 模 
色 填 充 算 法 处 理 样 本 中 的 异常 数据 ， 并 利用 健康 


型 中 的 未 知 参数 ; 其 次 ， 基 于 WGM(1，1) 模 型 ， 提出 动态 前 向 后 向 灰 


预测 过 程 预测 设备 的 剩余 有 效 寿命 ; 最 后 ， 通 过 实例 分 析 对 模型 进行 评价 和 验证 。 


效 解决 异常 数据 的 问题 。 


关键 词 : 分 段 隐 半 马 尔 可 夫 模 型 ; 期 望 最 大 化 自 适应 估计 算法 ; 动态 前 向 后 向 灰色 填充 算法 ; 寿命 预测 


中 图 分 类 号 : TP181 doi: 10.3969/j.issn.1001-3695.2018.01.0025 


结果 表明 ， 提 出 的 设备 健康 预测 方 


Equipment health prognosis based on dynamic EM-SHSMM under abnormal data 


Wu Jianfei, Liu Qinming' 
(Business School, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: Aiming at the problem of remaining useful life prognosis under abnormal data during equipment degradation, this 
paper developed a prognostic method based on dynamic expectation maximization (EM) -segmented hidden semi-Markov model 
(SHSMM) . First, based on the SHSMM model framework, it used the expectation maximization algorithm to estimate the 
unknown parameters of the model. Secondly, to process the anomaly data in the samples, it proposed a dynamic forward- 
backward gray-fill algorithm based on WGM (1, 1) , and it carried out the equipment health prognosis. Finally, it used a case 
study to evaluate the performance of the model. The results show that the proposed method could effectively solve the problem 
of abnormal data. 
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0 ”引言 测 数据 户 。Carr 和 Wang 考虑 了 多 源 监测 数据 下 的 剩余 寿命 预 
测 问 题 ， 基 于 拓展 卡尔 曼 滤波 技术 ， 提 出 了 一 种 近似 的 剩余 寿 

随 着 现代 化 生产 制造 技术 的 蓬勃 发 展 ， 许 多 机 械 设 备 变 得 。”” 命 预测 方法 ， 实 例 验 证 结果 表面 多 源 信息 融合 能 够 有 效 降低 剩 
越 来 越 复杂 。 这 些 机 械 设 备 由 于 运行 环境 、 负 载 等 因素 的 影响 ， 余 寿 命 预测 的 不 确定 性 四 , Kharoufeh 等 考虑 了 复杂 运行 环境 下 
其 性 能 及 健康 状态 将 会 发 生 不 可 避免 的 退化 ， 进 而 造成 设备 最 。 设备 的 可 靠 性 预测 问题 ， 提 出 了 非 时 齐 马 尔 可 夫 变 化 和 非 时 齐 
终 的 失效 。 一 旦 发 生 由 失效 引起 的 事故 ， 所 造成 的 人 员 财 产 损 。 半 马 尔 可 夫 变 化 两 种 方法 来 描述 随机 环境 的 变化 外 ,Dong 和 He 


康 状态 、 预 测 设备 的 剩余 使 用 寿命 ， 并 据 此 确定 对 设备 实施 维 方法， 搭建 了 一 个 集成 框架 多 传感器 收集 的 信 


失 甚 至 环境 破坏 往往 是 不 可 估量 的 。 因 此 ， 及 时 评估 设备 的 健 ”在 隐 马 尔 可 夫 (HMM) 基 础 上 提出 了 隐 半 马尔 可 夫 预 测 (HSMM) 


号 进行 融 口 ， YE 


护 的 最 佳 时 机 ， 对 于 切实 保障 复杂 设备 的 运行 安全 性 、 可 靠 性 ，” 著 提高 了 健康 预测 的 准确 性 54。Peng 和 Dong 又 对 隐 半 马尔 
与 经 济 性 具有 重要 的 意义 口 。 可 夫 模 型 进行 了 优化 改进 ， 通 过 引入 基于 役 龄 的 老化 因子 ， 从 

关于 健康 预测 的 研究 ， 引 起 了 国内 外 学 者 的 广泛 的 关注 ， 而 准确 描述 和 预测 生产 设备 的 性 能 劣化 中 。 针 对 设备 剩余 寿命 
并 在 最 近 十 年 得 到 了 长 足 的 发 展 。 例 如 Peng 和 Tseng 考虑 了 预测 问题 ，Liu 提出 了 一 种 基于 自 适 应 隐 式 半 马 尔 可 夫 模 型 
Wiener 过 程 的 漂移 系数 为 随机 变量 时 的 退化 建 模 吕 任 导 出 。 (AHSMM) 的 设备 健康 预测 方法 , 并 通过 实验 证 明了 该 方法 比 传 


题 ， 
了 寿命 分 布 的 显示 解 ， 但 该 方法 没有 利用 设备 运行 过 程 中 的 监 ” 统 的 HSMM 更 有 效 外 。 FangX 等 人 针对 退化 


言 号 缺失 问题 ， 提 
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全、 


则 研究 


出 了 一 个 半 参 数 方法 来 预测 部 分 退化 系统 的 剩余 寿命 由。 然而 ， ” 计 方 法 来 估计 6。 令 0={ 01, 02,…, 01 } 表 示 当 前 采集 到 的 观测 
前 的 剩余 寿命 预测 研究 大 都 假设 样本 数据 是 完整 的 0029。 但 值 ，5={ s,s2,…, ssi} 表示 观测 值 对 应 的 状态 ， 基 于 观测 值 0 的 


在 工程 实际 中 ， 由 于 噪声 、 扰 动 和 人 为 因素 的 影响 会 产生 不 合 。 对 数 似 然 函数 可 以 表示 为 
常理 的 数据 ， 而 在 设备 剩余 寿命 预测 研究 中 对 于 异常 数据 的 处 L(0)=1log[p(O/0)] (5) 
里 大 都 是 简单 剔除 或 者 不 处 理 (5-191。 其 中 :p(0/9) 为 观测 值 {07,02,…,o1} 的 联合 概率 密度 函数 ， 由 于 


异常 数据 的 存在 会 给 寿命 预测 带 来 较 大 的 偏差 ， 这 样 带 偏 。 5 是 一 个 隐 含 变量 , 因此 直接 使 上 式 关 于 9 最 大 化 是 不 可 行 的 。 
差 的 寿命 预测 应 用 于 健康 管理 时 势必 会 对 决策 结果 产生 影响 。 因此 ， 通 过 最 大 化 似 然 函 数 p(0,5/9) 并 对 隐 含 变量 8 求 取 关 于 
[1。 因 此， 针对 样本 数据 有 异常 的 情况 ， 本 文 在 SHSMM 的 整 ” 0,5/9 的 期 望 来 估计 逼近 参数 。 具 体 过 程 为 
体 架构 下 ， 首 先 利 用 期 望 最 大 化 自 适 应 估计 算法 估计 模型 的 未 E 步 又 : 计算 
知 参 数 。 其次, 在 改善 的 WGM(1,1) 模 型 基础 上 提出 动态 前 向 后 
向 灰色 填充 算法 对 样本 中 的 异常 数据 进行 处 理 ， 并 将 处 理 后 的 
完整 数据 输入 到 训练 好 的 SHSMM 模型 中 进行 状态 识别 与 寿命 。” 其 中 表示 在 第 i 步 估 计 的 参数 值 。 式 (2) 可 以 推导 为 


l 


Lr(0/ 0)=E S10,0 {log p(O.S /0)} (6) 


预测 。 L(916')=5p(0,S/ 6)ogp(0,510) 0) 
1 SHSMM 基本 理论 又 
结构 
SS P(LO.S10) = FE , (d= 4 qo)bs (OF )… 
与 常规 HMM 中 每 个 隐藏 状态 只 对 应 着 单一 的 观测 值 不 同 ， p (d= )» (0 ) (8) 
dss, Ndn-iJbs it 
SHSMM 中 每 个 隐藏 状态 都 对 应 着 一 个 片段 的 观测 值 ， 即 一 个 a 
观测 值 序列 。SHSMM 考虑 了 状态 的 持续 时 间 ， 所 以 当前 状态 Nl 
We A et ,A log p(O,S /0)=1logxs +》 logas 十 
向 下 一 状态 进行 转换 ， 不 仅 与 当前 所 处 的 状态 有 关 ， 还 与 当前 四 ba (9) 
所 处 状态 的 持续 时 间 有 关 。SHSMM 的 模型 表达 式 可 以 描述 为 a Oh ee 
G8-(x, 4, B,D)， 其 中 : x 为 初始 状态 概率 分 布 ，4 为 宏观 状态 转 下 
一 移 概 率 和 矩阵 ，B 为 观测 值 概率 矩阵 ， 忆 为 状态 驻 留 分 布 。 与 常 A 
T= 下 术 这 和 Ts nia 故 L(9/ 台 ) 可 以 表示 为 
中 ” 规 HMM 类 似 ，SHSMM 也 需要 解决 评估 、 识 别 和 训练 问题 。 
， ”1.2 推理 过 程 L(010)= Bogs, p(O0,S/0')+ 
首先 定义 一 个 前 向 变量 : 
&(i)=P(0.0,,…,0.,g=i/N) (1) 5 (Yow, 5S, ,Jo(0.s10)+ 
表示 在 t 时 刻 产 生 观 测序 列 (O1O2, …, On)， 并 终止 于 状态 i G ) 
logb, (0O% ,) |p(O,S/O')+.+ 
假设 设备 在 状态 i 持续 的 时 间 为 4， 概 率 为 P(dj)。0102， 


5 (>on, (d=4g, -4 ) jp(0.5/ 0) (10) 


n=1 


…, Or 表示 在 td 时 刻 以 状态 ;结束 而 产生 的 一 系列 观察 变量 
OraOrat1，…，0O1 表 示 在 1 时 刻 以 状态 j 结束 的 观察 变量 ， 将 在 M 步骤 : 计算 
f-d 时 刻 的 所 有 状态 和 这 些 状 态 的 可 能 持续 时 间 相 加 , 可 以 得 到 
以 下 递归 公式 : 


OD =argmax{L(0/0)) (11) 


WD 1 由 于 要 极 大 化 的 参数 单独 出 现在 方程 (6) 的 各 个 部 分 中 ,所 
wo(J)= wa_y(DwP(d) DO) 2 
和 l 1 | 以 只 需 对 各 个 部 分 分 别 极 大 化 。 利 用 拉 格 朗 日 乘 数 法 可 得 到 第 
其 中 :D 是 每 个 状态 可 能 持续 的 最 长 时 间 。 给 定 模型 4， 产生 观 “i+1 步 估计 的 参数 值 : 
察 序列 O 的 概率 为 


本 [| CDR(COTT CO) 


N 12) 
PolD=>wrO) G) DDS map aT OO 
与 前 向 变量 类 似 ， 后 向 变量 可 表示 为 jl 
os” 入 i MAGE HA) ACA A ER 
Bb(i)= 3p Wan Tl (0,) (4) 1 ww PII,(0.)B,()) 
2 ”期 望 最 大 化 参数 自 适 应 估计 算法 Bh nt (DB (D7) i 


为 了 解决 SHSMM 中 的 模型 训练 问题 ， 采 用 期 望 最 大 化 佑 并- 过-200401) 


录用 稿 


PAO AC IAN) YW RA 
DAOLI AAC NE) Y RA 


(15) 
然后 迭代 步骤 和 M 步骤 直到 p(O/B) 收 化 终止 ， 由 此 得 
到 对 应 的 参数 估计 值 。 


3 ”基于 SHSMM 的 设备 健康 预测 


3.1 WGM(1,1) 模 型 

在 传统 的 GM(1,1) 模 型 中 都 将 生成 系数 历 值 等 于 0.5， 然 
而 ， 这 并 没有 考虑 到 背景 影响 因素 对 预测 过 程 的 干扰 。 因 此 ， 
为 了 提高 模型 的 预测 精度 ， 对 原始 GM(1,1) 模 型 中 的 背景 值 进 
行 改进 ， 令 : 


P(d)= 


ZO(K)=WXO (K+ -WX KE-D) (16) 
将 ZO(k) 代 入 差分 方程 XX 中 (Kk)+aZ"(k)=b 中 ,并 通过 数 
归纳 法 可 得 预测 模型 为 : 
| 
(1+aW) 
利用 最 小 二 乘法 可 求 出 参数 a, bp, 再 将 a, bp, WW 的 值 代入 
到 方程 (16) 中 就 可 得 到 预测 值 XO(k) 。 
于 给 定 一 个 样本 数据 ,模型 的 精度 只 取决 于 生成 系数 到， 
所 以 如 何 找到 最 优生 成 系数 是 关键 。 令 &, =X®(K)-X(k)， 


3 


二 


,k=2,3... (17) 
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这 样 根 据 真实 值 与 预测 值 的 离 差 平方 和 5 = 六 妆 最 小 来 找到 使 


模型 精度 最 高 的 生成 系数 所 '， 并 最 终 以 该 生成 系数 下 的 
WGM(1,1) 模 型 进行 预测 。 
3.2 动态 前 向 后 向 灰色 填充 算法 

采集 设备 在 各 个 监测 时 刻 点 的 健康 状态 值 构成 样本 数据 自 
0=(XO(1)， 了 (2)，…， 了 0(n))。 动态 前 向 后 向 灰色 填充 算法 
的 基本 步骤 如 下 : 

a) 找 到 样本 数据 集中 的 第 一 个 异常 值 XO(k1) 并 将 其 标记 为 
缺失 值 ， 其 中 五 表示 第 1 个 异常 数据 的 监测 时 刻 点 。 

b) 在 COCOD，…，XoO(-D) 前 向 序列 基础 上 建立 WGMI(1,1) 
模型 ， 并 计算 相对 误差 ， 然 后 动态 的 向 前 增 大 r 个 步 长 ， 在 
COGc-D，…， XO(-1D)) 基 础 上 建立 新 的 WGM(1,1) 模 型 ， 并 计 
算 相 对 误差 。 循 环 上 述 步 又， 直到 满足 = 小 于 预 设 值 w。 选 择 
产生 相对 误差 最 小 的 前 向 灰色 预测 模型 ， 并 计算 出 异常 数据 前 
向 填补 值 六 (kh) 。 

co) 在 CCO(U+HD，…， oO)) 后 疝 序 列 基础 上 建立 WGM(1,1) 
模型 ， 并 计算 相对 误差 ， 然 后 动态 的 向 后 增 大 r 个 步 长 ， 在 
CO(+D，…，ZO( DJ) 基 础 上 建立 新 的 WGM(1,1) 模 型 ， 并 
计算 相对 误差 。 循 环 上 述 步骤， 直到 满足 rz 小 于 预 设 值 w。 选 
择 产 生 相 对 误差 最 小 的 后 向 灰色 预测 模型 ， 并 计算 出 异常 数据 


AT 


后 向 填补 值 x (4) 。 


d) 将 步骤 bjc) 中 的 计算 结果 进行 加 权 平 均 得 到 : 


ChinaX ivEé£ 
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要 ~ (0) 
RO)=— OK)+ EX (6) (18) 
忆 士 时 y+2z 


其 中 : y 为 最 优 前 向 灰色 模型 中 样本 数据 的 个 数 , z 为 最 优 后 向 


sy 


色 模 型 


补 值 。 


又 1 
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将 通过 


Pp 档 


9) 检 查 填 
否则 算法 终止 。 
3.3 ”健康 预测 过 程 


上 进行 设备 的 寿命 预测 。 


本 数据 的 个 数 ， 并 将 这 O(6) 作为 异常 数据 的 十 


后 的 样本 数据 是 否 存在 异常 值 ， 


若 存 在 返回 步 


动态 前 向 后 向 灰色 填充 算法 处 理 过 的 完整 观测 向 量 
输入 到 各 个 状态 下 的 SHSMM 中 ， 并 计算 观测 向 量 在 不 同 模型 
下 的 似 然 概率 对 数值 ， 其 中 输出 似 然 概率 对 数值 最 大 的 状态 即 
为 当前 设备 所 处 的 状态 ， 从 而 实现 设备 的 状态 识别 ， 在 此 基础 


和 和 


PE 


设备 的 全 寿命 历史 数据 训练 SHSMM 


Y 


获得 状态 转移 概率 和 各 个 状态 的 驻 留 时 间 
均值 和 方差 


当前 得 到 的 观测 值 


条 
六 
衣 

Ik 


和 


过 


动态 前 向 后 向 灰色 填充 算法 处 理 异常 
数据 
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识别 设备 当前 所 处 的 健康 状态 
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通过 方程 (19) 计 算 设 备 的 剩余 寿命 


图 1 健康 预测 基本 过 程 


基本 预测 过 程 如 图 1 所 示 。 


假设 设备 当前 所 处 的 健康 状态 为 状态 /，RUZ; 表示 设备 处 
于 健康 状态 7 的 平均 剩余 有 效 寿 俞 ， 则 
RUL = a [D(h)+ RUL, |]+ a (RUL,,) (19) 
案例 分 析 
为 了 验证 本 文 所 提 方 法 在 设备 状态 识别 与 寿命 预测 中 的 有 


效 性 


与 实 上 


用 性 ， 


数 所 


本 研究 应 用 美国 卡特 彼 勒 公司 液压 泵 的 全 寿命 
居 进 行 了 实验 分 析 。 液 压 泵 可 分 为 好 、 中 、 差 、 
在 实验 过 程 中 ， 针 对 每 一 种 条 件 收集 振动 信号 ， 每 个 状态 条 件 


坏 4 种 状态 。 


下 各 获取 30 个 数据 样本 ， 采 用 前 20 组 数据 进行 模型 训练 ， 后 


10 组 数 


4.1 


限定 在 e=0.000 1 。 
最 大 似 然 概率 的 对 数值 也 在 增加 ， 直 到 满足 收敛 误差 为 止 。 
于 SHSMM 的 训练 曲线 如 图 2 所 示 。 图 2 可 以 看 出 ，4 个 模型 


的 迭代 


不 超过 


线 均 可 以 达到 训练 时 设 定 的 误差 标 ; 
20 步 ， 另 外 在 利用 matlab 计算 过 程 中 ， 


有 于 测试 模型 ， 具 体 实验 数据 如 表 1 所 示 。 
模型 训练 
实验 中 设 定 SHSMM 训练 的 最 大 迭代 步 数 为 30, 收敛 误差 
在 模型 训练 过 程 中 ， 随 着 迭代 次 数 的 增加 ， 


[ou 


， 其 训练 步 数 均 


从 和 锡 代 开始 到 


收敛 结束 所 花费 时 间 不 超过 30s, 说 明 该 方法 训练 速度 快 、 训练 


精度 高 。 
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表 1 部 分 实验 数据 训练 好 的 模型 中 进行 测试 ， 测 试 结果 如 表 5 所 示 。 从 表 5 可 以 
状态 好 中 差 坏 看 出 异常 数据 会 降低 模型 的 识别 率 ， 而 处 理 后 的 SHSMM 识别 
1 2.62 2.41 31.26 41.78 率 为 100%。 
2 2.80 2.46 20.66 44.12 
3 2.65 2.58 8.89 31.08 表 4 SHSMM 状态 识别 结果 
4 3.22 2.37 5.06 57.04 健康 状态 好 中 差 坏 诊断 率 
5 2.39 2.53 41.07 60.15 好 9 1 0 0 90% 
6 2.38 2.42 7.53 108.03 二 0 这 
7 2.07 2.44 12.67 37.27 HMM 本 
8 2.24 2.54 5.74 35.61 关 0 1 8 1 80% 
9 2.25 2.58 29.22 69.00 坏 0 1 2 8 70% 
10 2.35 22.37 3.46 43.44 好 10 0 0 0 100% 
11 2.24 2.53 2.62 55.00 eo dd 005% 
12 2.28 2.46 3.66 44.12 SHSMM 
差 0 0 10 0 100% 
坏 0 0 0 10 100% 
了 表 5 异常 数据 下 和 处 理 后 的 状态 识别 结 
是 -| 健康 状态 好 中 差 坏 诊断 率 
E33 
起 好 5 1 0 0 833% 
Ei 中 0 4 2 0 667% 
异常 数据 下 的 SHSMM 
各 ao 差 0 1 3 50% 
KK 
坏 0 0 1 5 83.3% 
-250 5 10 有 20 25 30 好 6 0 0 0 100% 
基干 4 村 刑 站 | 结 中 0 6 0 0 100% 
图 2 ”基于 SHSMM 的 模型 训练 处 理 后 的 SHSMM 
差 0 0 6 0 100% 
利用 液压 泵 的 全 寿命 历史 数据 训练 得 到 各 个 状态 之 间 的 转 坏 00 0 6 100% 
移 概 率 以 及 每 个 状态 下 的 驻 留 时 间 均 值 和 方差 , 见 表 2 和 表 3。 
表 2 四 个 状态 之 间 的 转移 概率 43 大作 现 济 
状态 类 别 好 中 差 坏 通过 对 液压 泵 的 全 寿命 历史 数据 进行 训练 ， 可 以 获得 健康 
好 0 0s ee oe 状态 间 的 转移 概率 以 及 每 个 健康 状态 驻 留 时 间 的 均值 和 方差 。 
中 00d 人 和 Bo 且 经 过 测试 数据 检验 ， 通 过 训练 得 到 的 参数 值 都 是 可 行 的 。 
00000 5 0.9267 0.0733 方程 (19) 可 知 ， 设 备 的 剩余 有 效 寿命 等 于 当前 所 处 状态 下 驻 
。 v.0000 Go a L006 留 时 间 均 值 加 上 设备 下 一 个 健康 状态 下 的 剩余 有 效 寿命 。 根 据 
状态 识别 结果 确定 液压 泵 当前 所 处 的 健康 状态 ， 然 后 利用 方程 
& 3 ”四 个 状态 下 的 驻 留 时 间 均 信和 方差 (19) 计 算出 液压 泵 的 平均 剩余 有 效 寿 命 。 假 设 液压 泵 当前 处 于 
状态 类 别 好 中 关 坏 的 健康 状态 类 别 为 “ 差 ” 即 第 3 种 状态 。 由 表 2 和 表 3 知 ， 
均值 Osa ep a ee qa33=0.9267，D(h3)=11.4432，RUL4=10.5366, 计算 出 液压 泵 的 平 
设备 处 于 当前 状态 的 实际 剩余 寿命 21 的 相对 误差 仅 为 0.67% 1， 
4.2 状态 诊断 因此 ， 剩 余 预 测 结果 也 是 可 行 的 。 
将 各 个 状态 下 的 测试 数据 分 别 输入 到 训练 好 的 HMM 和 5 ”结束语 
a et ER Se 所 木 } 
SHSMM 模型 中 进行 测试 , 可 以 获得 40 个 测试 数据 的 状态 识别 i 
结果 ， 结 果 如 表 4 所 示 。 从 表 4 可 以 看 出 ， 基 于 SHSMM 的 状 本 文 提 出 了 一 种 针对 设备 健康 预测 中 异常 数据 如 何 处 理 的 
态 识别 率 为 100%， 明 显要 好 于 HMM 。 新 方法 。 提 出 的 动态 EM-SHSMM 模型 能 够 有 效 的 对 设备 退化 
分 别 从 各 个 状态 下 的 测试 数据 中 随机 选取 6 组 数据 ， 将 其 状态 进行 识别 以 及 剩余 寿命 的 预测 ， 采 用 期 望 最 大 化 自 适 应 估 
设置 为 异常 数据 ， 并 将 含有 异常 数据 的 测试 数据 输入 到 训练 好 计算 法 估计 模型 的 未 知 参数 , 基于 WGM(1,1) 模 型 , 提出 了 动态 


的 模型 中 进行 测试 ， 同 时 利用 提出 的 动态 前 向 后 向 灰色 填充 算 ”灰色 前 向 后 向 填充 算法 对 异常 数据 进行 处 理 ， 并 通过 对 卡特 彼 
法 对 异常 数据 进行 处 理 ， 然 后 将 处 理 后 的 完整 测试 数据 输入 到 ”” 勒 公司 液压 泵 的 状态 诊断 实验 证 明了 它 的 有 效 性 。 最 后 ， 根 据 
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本 文 所 提出 的 新 方法 能 够 有 效 解决 机 械 设 备 故 障 诊断 和 寿 
预测 中 存在 的 数据 异常 情况 。 但 是 , 本 文 仅 给 出 了 在 SHSMM 
的 设备 状态 诊断 和 寿命 预测 过 程 ， 而 健康 预测 的 目的 是 对 设 
进行 维修 决策 ， 因 此 接 下 来 的 研究 是 根据 SHSMM 的 预测 结 


吉 果 和 模型 的 训练 信息 可 以 得 到 设备 的 剩余 有 
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